GALEX Observations of CS and OH Emission in Comet 9P/Tempel 1 

During Deep Impact 1 

Paul D. Feldman 2 , Stephan R. McCandliss 2 , Jeffrey P. Morgenthaler 3 , 
Carey M. Lisse 4 , Harold A. Weaver 4 , and Michael F. A'Hearn 5 

pdf @pha . jhu . edu 
ABSTRACT 

GALEX observations of comet 9P/Tempel 1 using the near ultraviolet (NUV) ob- 
jective grism were made before, during and after the Deep Impact event that occurred 
on 2005 July 4 at 05:52:03 UT when a 370 kg NASA spacecraft was maneuvered 
into the path of the comet. The NUV channel provides usable spectral information in 
a bandpass covering 2000 - 3400 A with a point source spectral resolving power of 
R ~ 100. The primary spectral features in this range include solar continuum scattered 
from cometary dust and emissions from OH and CS molecular bands centered near 
3085 and 2575 A, respectively. In particular, we report the only cometary CS emis- 
sion detected during this event. The observations allow the evolution of these spectral 
features to be tracked over the period of the encounter. In general, the NUV emis- 
sions observed from Tempel 1 are much fainter than those that have been observed by 
GALEX from other comets. However, it is possible to derive production rates for the 
parent molecules of the species detected by GALEX in Tempel 1 and to determine the 
number of these molecules liberated by the impact. The derived quiescent production 
rates are g(H 2 0) = 6.4 x 10 27 molecules s _1 and Q(CS 2 ) = 6.7 x 10 24 molecules s _1 , 
while the impact produced an additional 1.6 x 10 32 H2O molecules and 1.3 x 10 29 CS2 
molecules, a similar ratio as in quiescent outgassing. 

Subject headings: comets: individual (9P/Tempel 1) — ultraviolet: solar system 
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1. Introduction 



On 2005 July 4 at 05:52:03 UT (as see n from Earth) a 370 kg NASA spacecraft was maneu- 
vered into the path of comet 9P/Tempel 1 (lA'Hearn et al.1 120051) . The results of the impact were 
recorded by a host of ground and space based o bservatories with wavelength coverage spanning 
the sub-mm to the ultraviolet (IMeech et al.ll2005l) . Here we report on the results of low resolution 
slitless spectroscopy acquired with the Galaxy Evolution Explorer (GALEX) in the near ultraviolet 
NUV channel 2 days before, immediately following, and 1 day after the impact event. The pri- 
mary spectral features in this range include a solar continuum scattered from cometary dust and 
emissions from OH and CS molecular bands centered near 3085 and 2575 A, respectively. 



2. Observations 

GALEX (Galaxy Evolution Explorer) is a NASA Small Explorer whose primary mission is 
to map the history of star formation using two modes: two-band photometry (FUV, 1350-1750 
A; NUV, 1750-3100 A) and integrated field grism spectroscopy with 10-20 A spectral resolution. 
Because of its large field-of-view, 1 .°2, GALEX is also well suited to cometary coma studies a s 



demonstrated by the 2005 March observations of C/2004 Q2 Machholz (jMorgenthaler et al hood) . 
It has a limiting spatial resolution in the NUV channel of « 5" sampled with 1"5 pixels. For Deep 
Impact only NUV data was obtained as the FUV side was not operating at time of the event. In all 
there were 7 contiguous pre-impact orbits devoted to grism observations, with « 900 s visibility 
per orbit. Two orbits of direct imaging were taken before and after the grism orbits. For the 
impact event there were 6 orbits of grism observations beginning a few minutes after impact, with 
direct imaging before (1 orbit) and after (2 orbits). One day post-impact there were 4 orbits of 
grism observations with one direct imaging orbit before the grism observations. A log of the 
spectroscopic observations is given in TableQ] Here we concentrate on the results derived from the 
low resolution slitless spectroscopy in the NUV channel. 

At the time of impact comet 9P/Tempel 1 was at a distance of r = 1.51 AU from the Sun and 
A = 0.89 AU from the Earth. A lf'5 pixel subtended 975 km at the comet. Emission from OH, CS, 
and solar scattered light were detected prior to impact. Just after impact the brightness and spatial 
extent of these emissions increased with time. The following day saw the return of these emissions 
to near the pre-impact levels. 
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2.1. Data Extraction 

The GALEX project supplied a time-tagged photon list of the observations with photon co- 
ordinates referenced to RA and DEC (J2000). The project also supplied the output of Emmanuel 
Bertin's SExtractor program run on the direct-mode (imaging) observations centered on the same 
fields. SExtractor provides accurate astrometry for sources found in the direct images. We con- 
verted these positions into rectangular strips using the GALEX dispersion relation for box length 
and the magnitude estimates from SExtractor for box widths in order to mask sources in the grism 
images. All of our astrometric corrections were done on a spherical projection with the IDLASTRO 
routines. Co-alignment between the strips and the grism images of the stars across the entire image 
confirmed the accuracy of these calculations. We used the JPL HORIZONS ephemeris generator 
with GALEX as the observatory ("@galex") to calculate the astrometrically correct ephemeris for 
comet Tempel 1. The ephemeris was used to reconstruct the comet image from the photon list 
in the reference frame of the comet. Counts and exposure times were accumulated into separate 
images and the source mask and flat field (supplied by the GALEX project) were applied every time 
the stars moved by one pixel relative to the comet. As a result, emission from well-characterized 
sources was effectively erased and very little exposure time was lost. The final rate images, con- 
structed by dividing the total count images by the total exposure time images for the corresponding 
observations, were used in the analysis of the temporal evolution of the coma following the impact. 

Because of the low total counts in each individual grism image, individual exposures were 
summed to produce a composite image from each day of the program. In this case the source 
mask was not applied and the stellar spectra appear as multiple images separated by the comet's 
motion between GALEX orbits. An example is shown in Figure \T\ for the sum of the six grism 
images obtained following the impact event. Note that the stellar images (each corresponding to 
a different orbit) are broadened by their apparent relative motion nearly normal to the dispersion 
direction during the exposure. 



3. Analysis 

Analysis of the spectrum is begun with the composite images (such as in Figure [Q) corre- 
sponding to the GALEX orbits from the three dates of observation. The spectra are extracted over 
13 rows of the rotated image array, corresponding to an effective slit height of 19. "5 to ensure that 
only pixels with statistically significant counts are included in the aperture photometry. The de- 
tector background is evaluated over the same rows in the spectral range 1700 to 2000 A where the 
NUV throughput for a solar-like spectrum is effectively zero. This background is subject to vari- 
ation due to the presence of faint stars or galaxies in the field. The three composite spectra, after 
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subtraction of background, are shown in the top panels of Figure [2] The extraction of the discrete 
emissions from CS and OH is complicated by the blended continuous solar spectrum produced by 
the scattering of sunlight by cometary dust grains in the coma. An additional complication arises 
from the different spatial distributions of the dust and gas components that result from the quies- 
cent outflow and the subsequent time evolution of these distributions after the impact, as the spatial 
variation is superimposed on the actual spectrum in the dispersion direction. 

Between 2000 and 3400 A the solar spectrum rises steeply and at the long wavelength end 
can contribute to the extended OH(0,0) band near 3085 A. Since the effective area, A e ff, supplied 
by the GALEX project did not extend beyond 3100 A, we identified a field F8 V star and an IUE 
calibration standard to extend A e ff to 3400 A. A model of the dust scattered solar continuum is 
required to account for the contamination of the OH(0,0) and C S(0,0) bands. We used a solar 
spectrum from UARS/Solstice, described by IWoods et al.1 (I1996Q . multiplied by this A e ff, and 
convolved with a superposition of ~80 A and ~250 A Gaussian instrument functions to simulate 
the radially outflowing dust in the dispersion direction, in order to generate a template for the 
solar scattered light. This template, normalized to re gions of the spectrum where strong molecular 
emission is absent, and reddened by 4% per 100 A (|Feldman & A'Hearnlll985r) . is also shown in 
the figure. The difference between the observed spectra and the template is shown in the lower 
panels of the same figure. 

The three identified emission features are the CS(0,0) band at 2576 A, and the OH(1,0) and 
(0,0) bands at 2820 and 3085 A, respectively. For near zero heliocentric velocity, the rat io of 
fluorescence efficiencies of the OH(1,0) and (0,0) bands is 0.10 (ISchleicher & A'Hearnlll988|) . but 
the GALEX effective area is 4.2 times higher at 2820 A than at 3085 A, so the (1 ,0) band should give 
a count rate ~0.4 that of the (0,0) band. This is consistent with the lower panels of Figure [2l and 
serves to validate the continuum subtraction. In principle, we could use the template scaling factor 
to derive the temporal variation of the dust ejected by the i mpact, but this is much better done with 



visibl e observations from the ground or from space (e.g., ISchleicher et al.l 120061 : Feldman et al. 
2007|) because of uncertainties in the background subtraction. 



GALEX observations began only a few minutes after the impact and grism spectra were 
recorded on six successive orbits, or for about 8.5 hours following the impact. The temporal 
evolution of the CS and OH emissions can be seen in the individual spectra shown in Figure [3] 
To construct a light curve, we used a rectangular aperture wide enough to capture almost all of 
the band emission, keeping the aperture height at 19. /; 5, as above. For CS, the width was 22."5, 
corresponding to a spectral bandwidth of 100 A, while for OH the width was 46 ."5, or ~200 A. 
The projected sizes on the sky are 12,600 x 14,500 km and 12,600 x 30,000 km, respectively 
for CS and OH. Similar photometry was performed on the pre- and post-impact data sets, and the 
results are shown in Figure @] in which the plotted times are the mid-points of the exposure. The 
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observations beginning about 24 hours after impact show the gaseous emissions returning to near 
quiescent levels. 



4. Discussion 



The quiescent production rates were determined using the vectorial model of iFestoul (|198ll) . 
The parent lifetimes at 1 AU were 71,000 s for H2O, with 85% of the dissociations producing OH 
(budzien et all 1994b . and 1 ,000 s for CS 2 (IFeldman et all 1999b . The lifetimes of OH and CS were 
taken to be 140,000 s and 100,000 s, respectively, but are not critical for the aperture size used. 
The parent outflow velocity in the model is 0.85r~ 1//2 km s -1 , while the daughters are given 1.05 
and 0.90 km s _1 , for OH and CS, respectively. The fluorescence e fficiencies at 1 AU are 2.3 x 
10 4 and 5.4 x 10 4 photons s 1 molecule 1 for the OH(0,0) band dSchleicher & AHearrJfl988b 
and the CS(0,0) band (R. Yelle, private communication), respectively. The derived production 
rates are given in Table |2l together with other recently reported values. Our value for the water 
production rate fits squarely within the range reported previously. We also note that the CS2/H2O 
production rate ratio, 9 .5 x 10~ 4 , is typical of values found from IUE observations of comets 
dMeier & AHearnlll997b . 



The linear increase with time following impact in the nu mber of dissociation p roducts in the 
aperture can be modeled, using known photochemical rates (IFeldman et al.l I2004Q . to determine 
th e abundance of the pa rent molecular species produced in the collision. This approach was used 
by iKtippers et al.l (12005b for the analysis of wide-field camera images of OH emission. However, 
for our aperture height of 12,600 km, molecules flowing radially outward at a mean velocity of 
0.5 km s -1 will begin to exit the aperture 3.5 hours after impact. This problem was considered 
by iManfroid et al.l (12007b . and we adopt their formalis m. In addition to the parent and daughter 
species, whose photochemical lifetimes are given above. IManfroid et al.l also include a "grandpar- 
ent" molecule to account for gas excavation over an extended period (several minutes) or evapora- 
tion from excavated icy grains. The "lifetime" of the grandparent is an adjustable parameter. 

The model results are shown with the data in Figure |U For a baseline, we adopt the mean 
value of pre-impact count rates, which correspond to an observation almost exactly one rotation 
period (~41 hours) earlier. The CS light curve is well fit by a mean velocity of 0.4 k m s _1 and 
a velocity dispersion of 0.1 km s~ ! , similar to the values derived by IManfroid et al.l for various 
radicals, and a "grandparent" lifetime of 2000 s. For OH, models with 0.4 km s _1 and 0.5 km s _1 
are shown in the figure. What matters most in deriving the number of parent molecules produced 
is the linear increase immediately follo wing the impa c t befo re any of the dissociation products 
has left the field-of-view. We note that IFeldman et al.l (I2006T) found a maximum velocity of CO 
produced by the impact to be 0.7 km s . 
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The total number of molecules found is ~ 1.3 x 10 29 molecules of CS2 and ~ 1.6 x 10 32 
molecules of H2O (corresponding to ~ 1.6 x 10 4 kg and ~ 4.6 x 10 6 kg, respectively). The 
CS2/H2O ratio is 8.0 x 10~ 4 , comparable to the quiescent value within the given uncertainties. 
A comparison with other reported values is given in Table[3j and shows generally good agreement. 
Our numbers are lower, by 25% for H2 O and a factor of 3 for CS2, than those presented in 2006 



in a preliminary report on these results (|Feldmanll2009l) . The spectral images used in that analysis 



did not have the comet's motion removed and so a large extraction box was used to ensure that all 
of the cometary photons were counted. This made it difficult to properly subtract the dust contin- 
uum, which is more important for isolating the CS emission than it is for the OH emission. We 
note that the only uncertainties in the derivation of the total number of molecules produced (other 
than due to photon statistics and absolute calibration) are those in the photodissociation lifetimes 
and fluorescence efficiencies, which are known at the 10% level. The combined uncertainties are 
^20%. 

After 25 hours all gas and dust emissions within the summation aperture have returned to 
nearly quiescent values. However, as the lower right panel of Fig. [2] shows, the OH bands one day 
after impact appear much broader in wavelength than prior to impact, likely due to outflow along 
the dispersion direction. This is not surprising as at 1.5 AU the photodissociation lifetimes of H2O 
and OH are 2 and 3 days, respectively. 

Two other orbiting spacecraft made NUV observations of the Deep Impact event: both are X- 
ray observatories with UV/optical imagers that were used with wide-band NUV filters. The filter 
passband of the Sw ift imager included both the CS and OH bands in addition to dust continuum 



(fMason et al J 120071) but the photometric images are dominated by the OH emission. The number 
of water molecules produced by the impact derived from their photometry (see Table [3]) is in 
excellent agreement with our result. The XMM-Newton observatory imager had two NUV filters 



one o f which included only CS emission while the other included both CS and OH (ISchulz et al. 



20061) . The light c urves derived from both filters show behavior similar to that of the two light 



curves in Figure |U ISchulz et al.1 attribute the difference in shape of the light curves to the presence 
of small "icy grains" in the ejecta. However, as the GALEX spectra of Figure [2] demonstrate, their 
data can be reasonably accounted for by CS emission. 



5. Conclusion 

GALEX observations of comet 9P/Tempel 1 at the time of Deep Impact provide a unique 
measurement of the CS parent produced by the excavation of cometary material, and as with other 
daughter molecules, appears to have the same abundance relative to H2O as in quiescent outgassing 
of the comet. Our derived value for the amount of H2O excavated is in excellent agreement with 
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the results of several other investigations. This reinforces the idea that nearly all the volatile species 
are uniformly mixed with the water ice from the near-surface layer at which sublimation normally 
occurs to the tens of meters to which Deep Impact excavated material. This lack of differentiation 
from near- surface to tens of meters implies either that the surface erodes (due to hydrodynamic 
drag on the dust by the gas) at about the same rate as the differentiation front proceeds deeper 
or that the differentiation has reached depths of many tens of meters while some how still leavin g 
the same relative fraction of all the volatiles throughout the upper tens of meters (lA'HearnlbOOSI) . 
Notwithstanding the relative faintness of the comet, the GALEX results demonstrate the unique 
capability of wide-field imaging spectroscopy of comets in the ultraviolet at wavelengths below 
the atmospheric cutoff. 



We thank the GALEX project, and Karl Forster in particular, for their efforts in planning 
and executing these observations. We would like to acknowledge Tim Conrow and the GALEX 
pipeline processing team for helping us to fully utilize the moving target capability of GALEX, Jon 
Giorgini and the JPL HORIZONS team for calculating high precision comet ephemerides from 
the perspective of GALEX, and Wayne Landsman and the many contributors to the IDLASTRO 
library. This work was supported by the GALEX Guest Investigator program under NASA grant 
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Table 1 . GALEX Grism Observation Log of Comet 9P/Tempel 1 . 



FILENAME EXPTIME (s) DATE-OBS TIME-OBS 



GI1_ 


_024001_ 


_DEEP_ 


_IMPACT_ 


_PRE_0001-ng-int. 


fits 


820.0 


2005- 


-07- 


-02 


06:22:17 


GI1_ 


_024001_ 


_DEEP_ 


_IMPACT_ 


_PRE_0002-ng-int. 


fits 


940.0 


2005- 


-07- 


-02 


08:00:53 


GI1_ 


_024001_ 


_DEEP_ 


_IMPACT_ 


_PRE_0003-ng-int. 


fits 


940.0 


2005- 


-07- 


-02 


09:39:29 


GI1_ 


_024001_ 


_DEEP_ 


_IMPACT_ 


_PRE_0004-ng-int. 


fits 


960.0 


2005- 


-07- 


-02 


11:18:04 


GI1_ 


_024001_ 


_DEEP_ 


_IMPACT_ 


_PRE_0005-ng-int. 


fits 


940.0 


2005- 


-07- 


-02 


12:56:41 


GI1_ 


_024001_ 


_DEEP_ 


_IMPACT_ 


_PRE_0006-ng-int. 


fits 


940.0 


2005- 


-07- 


-02 


14:35:17 


GI1_ 


_024001_ 


_DEEP_ 


_IMPACT_ 


_PRE_0007-ng-int. 


fits 


940.0 


2005- 


-07 


-02 


16:13:53 


GI1_ 


_024002_ 


_DEEP_ 


_IMPACT_ 


.0001-ng- int. fits 




759.0 


2005- 


-07 


-04 


06:01:47 


GI1_ 


_024002_ 


_DEEP_ 


_IMPACT_ 


_0002-ng-int.fits 




959.0 


2005- 


-07 


-04 


07:40:23 


GI1_ 


_024002_ 


_DEEP_ 


_IMPACT_ 


_0003-ng-int.fits 




939.0 


2005- 


-07- 


-04 


09:18:59 


GI1_ 


_024002_ 


_DEEP_ 


_IMPACT_ 


_0004-ng-int.fits 




938.1 


2005- 


-07- 


-04 


10:57:36 


GI1_ 


_024002_ 


_DEEP_ 


_IMPACT_ 


_0005-ng-int.fits 




940.0 


2005- 


-07- 


-04 


12:36:11 


GI1_ 


_024002_ 


_DEEP_ 


_IMPACT_ 


_0006-ng-int.fits 




939.0 


2005- 


-07 


-04 


14:14:48 


GI1_ 


_024003_ 


_DEEP_ 


_IMPACT_ 


_POST_0001-ng-int 


.fits 


299.0 


2005- 


-07 


-05 


05:05:56 


GI1_ 


_024003_ 


_DEEP_ 


_IMPACT_ 


_POST_0002-ng-int 


.fits 


919.0 


2005- 


■07- 


-05 


06:40:52 


GI1_ 


_024003_ 


_DEEP_ 


_IMPACT_ 


_POST_0003-ng-int 


.fits 


931.2 


2005- 


-07 


-05 


08:19:29 


GI1_ 


_024003_ 


_DEEP_ 


_IMPACT_ 


_POST_0004-ng-int 


.fits 


939.0 


2005- 


-07- 


-05 


09:58:05 



Table 2: Quiescent production rates (in molecules s ) in 9P/Tempel 1. 



Reference 


Observation 


H 2 O(xl0 27 ) 


CS 2 (xl0 24 ) 


This paper 


OH/CS - UV (GALEX) 


6.4 


6.7 


Keller et al. (2001) 


OH - UV (Rosetta) 


4.2-4.4 




Schleicher et al. (2006) 


OH - UV (ground) 


6 




DiSanti et al. (2007) 


H 2 - IR 


9.4 




Bensch et al. (2007) 


H 2 - sub-mm (SWAS) 


5.2 




Biver et al. (2007) 


H 2 0/CS - sub-mm (ODIN) 


9.1 


< 19.7 a 



"Post-impact on 6 July. 



Table 3: Total number of water molecules produced by the impact. 



Reference 


Observation 


H 2 O(xl0 32 ) 


This paper 


OH 


- UV (GALEX) 


1.6 


Mason et al. (2001) 


OH 


- UV (SWIFT) 


1.4±0.2 


Keller et al. (2007) 


OH 


- UV (Rosetta) 


1.5-3.0 


Bensch et al. (2007) 


H 2 0- 


sub-mm (SWAS) 


< 5.9 (3-G) 


Biver et al. (2007) 


H 2 0- 


sub-mm (ODIN) 


1.7±0.7 
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FIGURE CAPTIONS 

Fig. 1. — The central 12/8 x 12/8 of the composite GALEX NUV grism image taken in the 8.5 h 
following the Deep Impact encounter with comet 9P/Tempel 1 . The motion of the comet is nearly 
perpendicular to the dispersion as can be seen from the multiple, wide stellar spectra. The image 
has been rotated 48° east from north. Nearly point-like CS and continuum emission together with 
diffuse OH emission at 3085 A can be seen in the spectrum. 

Fig. 2. — Composite GALEX NUV spectra of comet 9P/Tempel 1 for each of the three days of 
observation as listed in Table [TJ Background counts have been subtracted. The top panels show 
the fit to solar scattered light in red, while the lower panels show the spectrum with the solar light 
subtracted. The peak wavelengths of the CS (0,0) and OH (1,0) and (0,0) bands are indicated. 
Note, in the post-impact spectra, the broadening of the OH emissions due to outflow along the 
dispersion direction. 

Fig. 3. — Same as the top panels in Figure |2] for the six consecutive GALEX orbits taken immedi- 
ately following the impact, showing the temporal evolution of the emissions. 

Fig. 4. — Temporal evolution of the CS (top) and OH (bottom) emissions. The left panels show the 
count rates roughly one comet rotation (~41 hours) prior to the impact. Model fits to the aperture 
photometry are shown for velocities of 0.4 km s _1 (dashed line) and 0.5 km s _1 (dotted line). 
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Fig. 1. — The central 12/8 x 12/8 of the composite GALEX NUV grism image taken in the 8.5 h 
following the Deep Impact encounter with comet 9P/Tempel 1 . The motion of the comet is nearly 
perpendicular to the dispersion as can be seen from the multiple, wide stellar spectra. The image 
has been rotated 48° east from north. Nearly point-like CS and continuum emission together with 
diffuse OH emission at 3085 A can be seen in the spectrum. 
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Fig. 2. — Composite GALEX NUV spectra of comet 9P/Tempel 1 for each of the three days of 
observation as listed in Table [U Background counts have been subtracted. The top panels show 
the fit to solar scattered light in red, while the lower panels show the spectrum with the solar light 
subtracted. The peak wavelengths of the CS (0,0) and OH (1,0) and (0,0) bands are indicated. 
Note, in the post-impact spectra, the broadening of the OH emissions due to outflow along the 
dispersion direction. 
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Fig. 3. — Same as the top panels in Figure [2] for the six consecutive GALEX orbits taken immedi- 
ately following the impact, showing the temporal evolution of the emissions. 
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Fig. 4. — Temporal evolution of the CS (top) and OH (bottom) emissions. The left panels show the 
count rates roughly one comet rotation (~41 hours) prior to the impact. Model fits to the aperture 
photometry are shown for velocities of 0.4 km s _1 (dashed line) and 0.5 km s _1 (dotted line). 



